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ABSTRACT 

We present a new observational study of the 12 CO(1-0) line emission as an 
H2 gas mass tracer under extreme conditions in extragalactic environments. Our 
approach is to study the full neutral interstellar medium (H2, H I and dust) of two 
galaxies whose bulk interstellar medium (ISM) resides in environments that mark 
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(and bracket) the excitation extremes of the ISM conditions found in infrared 
luminous galaxies, the starburst NGC 3310 and the quiescent spiral NGC 157. 
Our study maintains a robust statistical notion of the so-called X = N(B. 2 )/Ico 
factor (i.e. a large ensemble of clouds is involved) while exploring its dependency 
on the very different average ISM conditions prevailing within these two systems. 
These are constrained by fully-sampled 12 CO(3-2) and 12 CO(1-0) observations, 
at a matched beam resolution of Half Power Beam Width ~ 15", obtained with 
the James Clerk Maxwell Telescope (JCMT) on Mauna Kea (Hawaii) and the 
45-m telescope of the Nobeyama Radio Observatory (NRO) in Japan, combined 
with sensitive 850 /iin and 450 /im dust emission and H I interferometric images 
which allow a complete view of all the neutral ISM components. Complementary 
12 CO (2-1) observations were obtained with JCMT towards the center of the two 
galaxies. We found an X factor varying by a factor of 5 within the spiral galaxy 
NGC 157 and about 2 times lower than the Galactic value in NGC 3310. In 
addition, the dust emission spectrum in NGC 3310 shows a pronounced submm 
"excess". We tried to fit this excess by a cold dust component but very low 
temperatures were required (T c ~ 5 — 11 K) with a correspondingly low gas- 
to-dust mass ratio of ~ 5 — 43. We furthermore show that it is not possible 
to maintain the large quantities of dust required at these low temperatures in 
this starburst galaxy. Instead, we conclude that the dust properties need to be 
different from Galactic dust in order to fit the submm "excess". We show that 
the dust SED can be fitted by an enhanced abundance of Very Small Grains 
(VSGs) and discuss different alternatives. 



1. Introduction 

The use of the 12 CO(1-0) line emission as a tracer of H 2 gas has become an established 
technique since its the first detection (Wilson, Jefferts & Penzias 1970), with several large 
surveys using it to obtain images of its distribution in galaxies across the Hubble sequence 
(e.g. Young & Scoville 1991; Heifer et al. 2003; Leroy et al. 2005) and examine its relation 
to star formation (e.g. Wong & Blitz 2002). Several theoretical studies have been published 
since in order to outline the dependence of the so called X factor converting the luminosity 
of this optically thick tracer to H 2 mass under the ambient ISM conditions (e.g. Dickman, 
Snell, & Schloerb 1986; Wolfire, Hollenbach, & Tielens 1993; Bryant k Scoville 1996). These 
have now converged towards a view of a fairly robust X factor provided that: a) the 12 CO(l- 
0) line is well-excited, b) a large ensemble of CO-bright molecular clouds is involved (so the 
statistical notion of X remains valid), c) the molecular gas observed is reducible to a set of 
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self-gravitating "units" (whose linewidths reflect solely their gas mass), and d) for a given 
far-ultraviolet (FUV) interstellar radiation field the molecular gas is sufficiently metal-rich 
to ensure that CO remains abundant throughout the volume of a typical molecular cloud 
(Z > 1 for G = 1 in Habbing units). The latter restriction stems from the fact that unlike 
H2 the CO molecule, with its much lower relative abundance ([CO/H2] ~ 10~ 4 ), cannot 
self-shield from the dissociating FUV radiation and its relative extent in molecular clouds is 
regulated by the absorption of that radiation by dust (Pak et al. 1998). 

The strong dependence of the X factor on metallicity has been expected from theo- 
retical studies (Maloney & Black 1988; Bolatto, Jackson, & Ingalls 1999), and confirmed 
by extensive observations of the Magellanic Clouds (Israel 1997, 1999) and the Magellanic 
irregular IC 10 (Madden et al. 1997). In the high-excitation ISM environments found in 
the Ultra Luminous Infrared Galaxies (ULIRGs) an unbound molecular gas phase (a result 
of their violent, merger-induced, kinematic environments, e.g. Aalto et al. 1995) is thought 
as responsible for a systematic overestimate of the H 2 gas mass when a Galactic X factor 
Xg is used (Solomon et al. 1997; Downes & Solomon 1998). Towards the other excitation 
extreme, a cold and/or diffuse gas phase with very low 12 CO(1-0) brightness has been dis- 
covered in molecular clouds in the nearby quiescent galaxy M 31 (Allen & Lequeux 1993; 
Allen et al. 1995; Loinard, Allen, & Lequeux 1995), as well as in the outer Galaxy, leading 
even to suggestions that cold molecular gas constitutes an important component of the dark 
matter in spirals (Lequeux, Allen & Guilloteau 1993). 

Unlike the relatively well-explored metallicity dependence of the X factor, the effects of a 
highly excited, possibly non self-gravitating, molecular gas phase in starburst environments, 
and a low-excitation one in quiescent ones harboring low star formation rates, have not 
been widely studied. For ULIRGs such efforts have been hindered by their compact CO 
emission regions, the presence of a "standard" ensemble of star-forming molecular clouds 
mixed with the highly-excited phase, and the current inability to conduct multi-line CO 
observations at high and beam-matched resolutions. The latter could disentangle the two 
gas phases present in those systems, and such observations of ULIRGs will be important in 
the future. In low excitation environments, identified by the presence of cold dust and/or low 
12 CO(2-1)/ 12 CO(1-0) line ratios (r 2 i), lack of data constraining the gas excitation allows a 
degeneracy between a diffuse and a cold gas phase, which translates to large uncertainties 
in the estimated H 2 gas mass (e.g. Israel, Tilanus, & Baas 1998). Moreover, while it can 
be argued that very cold molecular clouds do exist in a few places in galaxies, it is far from 
certain that they can be typical for a significant fraction of their total molecular gas mass, 
even in quiescent environments (e.g. Wilson & Mauersberger 1994; Sodroski et al. 1994). 

From the largest CO line ratio survey available (Braine & Combes 1992, BC92 hereafter) 
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we chose two infrared luminous galaxies, NGC3310 (a starburst) and NGC 157 (a quiescent 
spiral), whose physical conditions of the bulk of their molecular gas and dust mass best 
represent the excitation extremes within the IRAS Intermediate Luminosity Sample [Lm = 
(10 10 — 10 11 )L Q ] (Tinney et al. 1990). Both have similar far-infrared (FIR) luminosities 
(~ 10 10 L Q ) but very different S 100 /S 60 "colors" of ~ 1.28 (NGC 3310) and ~ 2.42 (NGC 157) 
while their reported 12 CO(2-1)/ 12 CO(1-0) line ratios mirror this with r 21 (NGC 3310)~ 2.6 
(the highest reported such ratio, but see §3.1.1), and r 2 i(NGC157)~ 0.4 (typical for cold 
quiescent molecular clouds in the outer Galaxy). We make a comparison study of their ISM 
properties using the 12 CO(3-2) and 12 CO(1-0) lines, along with 850 yum and 450 /xm dust 
continuum maps and sensitive H I images obtained from the literature. Our main goal is to 
investigate excitation effects on the H 2 gas mass estimates in the extragalactic domain. 

The optically thin dust continuum emission at FIR and submillimetre wavelengths can 
be used to measure the dust column density, and thus provide an alternative way to trace the 
gas mass in most ISM environments where gas and dust are well mixed. Such an approach 
is independent of gas excitation and can potentially detect very cold gas in regions where 
CO is not excited (Bernard et al. 2008), but its use in extreme extragalactic environment 
is hindered by a large uncertainty in the gas-to-dust (GD) mass ratio. We compute the GD 
mass ratios in our sample galaxies and explore whether we can assume a Galactic GD ratio 
to estimate the gas mass from dust continuum emission. 

In §2 we describe our observations and data reduction procedures. In §3 we present our 
results from the molecular line analysis as well as the dust SED modeling, while in §4 we 
discuss the extreme ISM environments that the two galaxies under study represent. Finally, 
in §5 we summarize our conclusions. Appendix A shows how the interstellar radiation field 
for the two galaxies was calculated. 

Throughout the paper we assume a Hubble constant of H Q = 50 km s _1 Mpc -1 . The 
distances in our calculations are assumed to be 21.3 Mpc for NGC 3310 and 35.0 Mpc for 
NGC 157 (Stevens et al., 2005). 
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2. 



Observations and data reduction 




We mapped the 12 CO (1-0) emission of both galaxies with the Nobeyama 45-m telescopqj 
and the 12 CO(3-2) emission with the James Clerk Maxwell Telescope (JCMT)@. Given the 
importance of accurate measurements of the 12 CO(3-2)/ 12 CO(l-0) line ratio as a function 
of position in our present study, we devoted a large amount of our observing time at both 
telescopes towards pointing checks and observations of spectral line standards. We also 
used the two telescopes to observe 12 CO(3-2) and 12 CO(1-0) towards a common position 
in Orion-KL. The physical conditions of the H2 gas towards this region are constrained 
by a multitude of molecular line observations reaching as high as 13 CO(5-4), with inferred 
densities of n(H 2 ) > 10 4 cm -3 and T^ in ~ 45 K (Plume et al. 2000). For such conditions the 
intrinsic 12 CO(3-2) / 12 CO(1-0) brightness temperature ratio expected is R31 ~ 1 and thus for 
the similar beams of the JCMT and the Nobeyama 45-m telescope: T m b(3 — 2)/T m b(l — 0) ~ 
i? 3 i ~ 1. Our calibration observations found very good agreement for associated spectra 
from both telescopes, with an averaged integrated intensity ratio of 0.93. This ensures that 
there are no serious calibration offsets between the two telescopes. 



We used the JCMT during the periods of July 1999, January 2000 and February 2001 to 
conduct 12 CO(3-2) imaging of the starburst galaxy NGC3310 and then again in June 2002, 
and May, June and July of 2004 for similar observations of the quiescent spiral NGC 157. 
In all cases we employed receiver B3 in a dual channel mode, tuned to Single Side Band 
(SSB) to observe the 12 CO(3-2) line at 345.796 GHz and used the Digital Autocorrelation 
Spectrometer (DAS) with a bandwidth of 920 MHz (~ 795 km s _1 ). Both maps were acquired 
with a A6 = 1/2 6'hpbw = 7" sampling, which resulted in the acquisition of spectra in 47 
positions for NGC 3310, and 194 positions in NGC 157, covering most of the dust emission 
detected in our submillimeter (submm) continuum maps of these galaxies. We employed 
rapid beam- switching at a frequency of 1 Hz and chop throws of ~ 120" — 180" (in Az). For 



1 The Nobeyama Radio Observatory (NRO) is a division of the National Astronomical Observatory of 
Japan (NAOJ) under the National Institutes of Natural Sciences (NINS). 

2 The JCMT is operated by the Joint Astronomy Center in Hilo, Hawaii, on behalf of the parent organi- 
zations Science and Technology Facilities Council in the UK, the National Research Council of Canada, and 
the Netherlands Organization for Scientific Research. 



2.1. 



CO observations 



2.1.1. 12 CO (3-2) and 



CO(2-l ) observations at the JCMT 
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the NGC3310 observations typical system temperatures where T sys ~ (400 — 600) K with 
At int ~ 20 min of integration time per point (effective time for dual channel operation), 
while for NGC 157 these figures were T sys ~ (300 - 400) K, and At int ~ (10 - 20) min. The 
pointing was monitored every (45-60) min and from the pointing residuals we estimated an 
rms pointing uncertainty of S9 Tms ~ 2" — 3". 

A position with strong 12 CO(3-2) emission was repeatedly observed within each galaxy 
during our imaging observations to verify proper pointing and consistent line calibration (i.e. 
by checking line profiles and intensities). These observations along with numerous ones of 
spectral line standards yielded a spectral line calibration uncertainty of ~ 10%. The latter 
is the 2nd term contributing to the total uncertainty of line intensities given by 




where T is the main beam brightness temperature of the line, averaged over a total number 
of channels N c ^, and with a baseline defined over a total number of channels N^. The 
second term describes the uncertainties inherent in the calibration of mm/submm lines (e.g. 
temperature uncertainties of cold/warm loads used in the 3-point calibration) and the 3rd 
term is the uncertainty associated with the adopted value of the main beam efficiency rj m ^. 
The contribution of the pure thermal noise is given by 



5T\ 5T ch fN has + N ch \ 1/2 



TVther T V N has N ch 



where 5T& is the thermal rms noise across the line-free channels of the spectral line. Repeated 
observations of Mars and Uranus yielded beam efficiency estimates at 345 GHz of 7/ m b = 0.62 
for the NGC 3310, and r/ m b = 0.65 for the NGC 157 observations. Inspection of the JCMT 
databaj demonstrates a ~ 10% uncertainty for the rj m b values. The latter enters the total 
rms uncertainty (equation 1) of the spectral line strengths in our maps (and thus that of the 
12 CO(3-2)/ 12 CO(l-0) (r 31 ) line ratios), but is not expected to contribute to the expected 
rms line strength (and thus line ratio) variations within maps. 

We also used receiver A3 tuned to Double Side Band (DSB) to observe the 12 CO(2- 
1) line (230.538 GHz) towards the center of both galaxies. The DAS was used at a 1.8 
GHz bandwidth (~ 2340 km s _1 ). Observations of Mars, Uranus and Jupiter yielded beam 



3 http:/ /www .jach.hawaii.edu/JCMT/spectralJine/Standards/efF_wcb. html 
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efficiencies of r] m b = 0.67 for NGC3310 and r) m b = 0.71 for NGC 157, for a beam of #hpbw ~ 
21", values that we adopt for our conversion of the T£ to the T mb scale at 230 GHz. 

Finally all the spectra per position were inspected individually and we co-added only 
those found to be consistent within the S/N and after removal of bad channels on a few 
occasions. Linear baselines were subtracted and the 12 CO(3-2) spectra were binned to a 
resolution of 25 MHz (~ 21.6 km s" 1 ), divided by the appropriate rj m b and then added to the 
grid maps of the galaxies (Figures 1 and 2). The maps were convolved to a final resolution of 
15" (as to exactly match that of the 45-m) and then inserted into AIPS for further processing 
and combination with the 12 CO(1-0) cubes. 

2.1.2. 12 CO(1-0) observations at the Nobeyama 45-m telescope 

The observations with the 45-m telescope of the National Radio Observatory (NRO) 
of Japan were made in December 2000 and December 2004. We used the Super conductor- 
Insulator-Superconductor (SIS) receiver S100 as the frontend and an Acousto-Optic Spec- 
trometer (AOS), with a frequency resolution of 250 kHz and the total bandwidth of 250 MHz 
(~ 650 km s _1 ). The intensity calibration was performed by the two-point chopper-wheel 
method, yielding the total T sys , which typically was T sys ~ (500 — 700) K. The telescope 
pointing was checked every 30 min by observing SiO maser emission of late type stars at 
43 GHz. The rms pointing uncertainty was found to be 59 rms ~ 4". We used a position- 
switching mode with an integration time of 20 sec for on- and off-source. We measured 
the main beam efficiency using planets and found it to be i] m b ~ 0.4 with an uncertainty 
of ~ 15%. The half-power beam width was estimated from observations of 3C 279 and is 
#hpbw ~ 15". We finally obtained several spectra of 12 CO(1-0) of Orion-KL in order to 
obtain the line calibration intensity as well as check the relative calibration of the Nobeyama 
45-m telescope with respect to the JCMT. The integration time was 10 min for most points, 
thus the rms is higher in the 12 CO(1-0) spectra, and the baseline is also worse than that of 
the 12 CO(3-2) due to position switching, especially for the weak 12 CO(1-0) lines. 

2.2. Dust continuum submm observations with SCUBA 

Our SCUBA observations of NGC 3310 and NGC 157 were carried out on the nights 2000 
January 24 and 25, with grade 1 and 2 weather respectively. The SCUBA submm bolometer 
array (Holland et al. 1999) is an array of 37 bolometers at 850 /xm and 91 at 450 /im covering 
a region of sky of ~ 2.3arcmin in diameter. In order to provide fully sampled images, the 
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secondary mirror moves in a 64-step jiggle pattern, with the integration time lasting Is at 
each position. After the 16 steps of jiggle pattern, the telescope nods in order to allow for 
slowly varying sky gradients. 

Additionally to the observations obtained by us, we used SCUBA data from the JCMT 
archive (see §3.3.1). The data were reduced using the SURF software package (Jenness et 
al. 1998). We made corrections for the atmospheric absorption using the opacities derived 
from skydip measurements taken regularly every night. Noisy bolometers were flagged and 
large spikes were removed. During each night, several calibration sources were mapped in 
order to flux calibrate our images. Finally, the individual jiggle maps were co-added, after 
the proper reduction was made, and a final map was created for each galaxy. From the 
450 /iin observations described above, only the ones taken in grade 1 weather conditions 
had sufficient signal-to-noise ratio and were used in our analysis. Flux calibration of the 
final maps, in units of Jy beam -1 , as well as the corresponding beam profiles were obtained 
from beam maps of CRL 618 observed with the same chop throw as our galaxy maps. The 
calibration information is summarized in Table 1. 

3. Results 

3.1. The CO data 

In Tables 2 and 3 we tabulate the integrated intensity of 12 CO(1-0), 12 CO(3-2) (7 10 , 
^32) and the r 31 = /32//10 ratios for the different positions within NGC3310 and NGC 157 
respectively. The uncertainty in iio, ^32 is calculated according to equation (1) and the 
uncertainty for r 3i was combined quadratically from that of iio and 1^2- In most cases a w 
is higher than <7 32 for the reasons described in §2.1.2, hence the uncertainty in r 31 is mostly 
from o" 10 . A pointing error could cause extra uncertainty in r 31 when the two telescopes were 
not observing the same position. The expected error is on average 10-20% for a pointing 
error of 4" at 12 CO(1-0) and 3" at 12 CO(3-2), but could be in principle as high as ~ 50% if 
both telescopes are mis-pointed in different directions by a full grid cell size of 7". 

In Figures 1 and 2 we present the 12 CO(3-2) (red) and 12 CO(1-0) (green) profiles 
overlaid on a grey scale SCUBA image of the dust emission at 850 /im of NGC 3310 and 
NGC 157 respectively. In all positions where high S/N 12 CO(3-2) and 12 CO(1-0) spectra 
exist there is good agreement in their line profiles, demonstrating a good relative pointing 
on the grid positions by the two telescopes. 
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3.1.1. NGC3310: a high excitation case 

For NGC3310, 12 CO(1-0) and 12 CO(3-2) emission is detected in 33 positions covering 
a 42" x 42" region, corresponding to 4.3 kpc x 4.3 kpc in physical scale (see Figure 1 and 
Table 2). The r 31 ratio is slightly higher than unity in most positions, and slightly enhanced 
mostly in the southwest in comparison to the northeast of the nucleus (though only marginal 
considering the S/N). The area average is (r 3 i) = 1.17, while the range of values found within 
the inner region of NGC 3310 is r 31 = 0.9 — 1.4. This range is still consistent with a constant 
ratio of (r 3 i) =1.15 ± 0.25 (the average value of the extreme r 3i values mentioned above) 
given the ~ 20% noise variations of this ratio expected within the map. These are very 
high values, suggesting a highly excited molecular gas phase throughout NGC 3310. For 
comparison, starburst nuclei show an average ratio of r 3 i = 0.64 (Devereux et al. 1994; Yao 
et al. 2003; Narayanan et al. 2005). 

In the survey conducted by Braine and Combes (1992) NGC 3310 had the highest r 21 
reported with r 21 = 2.6. This value is high even when compared to the highest r 21 values 
found within individual molecular clouds (molecular line ratios of individual regions within 
molecular clouds span a much larger range than such ratios over large molecular cloud 
ensemble averages encompassed by most extragalactic observations). Interestingly r 2 i > 
2 ratios have only been observed in the N83/N84 region of the Small Magellanic Cloud 
(Bolatto et al. 2003), which is part of an expanding molecular shell centered on a supernova 
remnant (SNR; Haberl et al 2000). For much larger scales such high ratios have been 
reported also for the nearby archetypical starburst galaxy M 82 (Loiseau et al. 1990), but 
were not confirmed by subsequent observations (Mao et al. 2000 and references therein). 
This underscores the need for beam-matched observations (that become progressively more 
difficult with the higher J-separation of rotational transitions) and good cross-calibration of 
mm/submm telescopes for accurate line ratio measurements. 

Our own 12 CO (2-1) measurements (Figure 3a) yield velocity-integrated line intensity 
I = J T-xfodV of I21 = 9.8 ± 1.6 K km s _1 for NGC 3310, in excellent agreement with the 
value I 2 \ = 9.5 Kkms" 1 reported by BC92. However for 12 CO(1-0) at the (0, 0) position, 
after convolving to a 21" resolution we find I w = 6.62 ± 2.00 Kkms -1 , which is 1.8 times 
higher than that reported by BC92. As a result, r 2 i is 1.48 ± 0.25 at the central 21" region, 
which is significantly lower than the value reported in BC92, but still indicative of the 
highly excited molecular gas that seems to pervade this galaxy. Convolving the 12 CO(3-2) 
map to a 21" resolution yields J 32 = 9.8 ± 1.8 and r 3 i = 1.48 for the (0",0") position. The 
extreme ISM excitation conditions in NGC 3310 revealed by the 12 CO(1-0), 12 CO(2-l), and 
12 CO(3-2) relative line strengths are further corroborated by the warm IRAS "color" ratio of 
Sioq/Sqq ~ 1.28 (compared, e.g. with the value of ~ 2.42 for the quiescent spiral NGC 157). 
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3.1.2. NGC 157: variations in the gas excitation 

For NGC 157, 12 CO(1-0) and 12 CO(3-2) emission is detected in 64 positions within a 
28" x 126" region, corresponding to 4.7 kpc x 21.4 kpc in physical scale (see Figure 2 and 
Table 3). The r 31 ratios are slightly higher than unity at the nuclear region but drop to less 
than 0.5 in the two spiral arms, in a remarkably different fashion from the r 31 distribution 
seen in NGC 3310. Indeed, quite unlike the starburst NGC 3310 whose galaxy-wide star- 
forming activity leaves no discernible gas excitation gradient, NGC 157 is quiescent, with 
little star formation activity but with a marked change of ISM physical conditions. Maps of 
12 CO(3-2) line emission in 12 nearby galaxies studied by Dumke et al. (2001) have shown 
such decreasing r 31 ratios outward from the centers of many galaxies but none are as large 
as those measured in NGC 157. 

High r 3 i ratios are found near the nucleus at offset (0",0") and (0", 7"), coincident 
with high SCUBA S^o/Ssm flux ratios, suggesting a high dust and gas temperature at this 
position. Low values of r 31 ~ 0.2 on the other hand are found in the northeast (NE) and 
southwest (SW) disk, near offset (7", 49") and (-7", -42") respectively, which host strong 
850 fim emission with a low S^o/Ssm flux ratio and thus low dust temperature (see §3.3.3). 

For the central 23" of this galaxy Braine et al. (1993) report iio = 23 ± 0.9Kkms _1 . 
However, after convolving our 12 CO(1-0) data to 23", we find iio = 10.7 ± 1.5Kkms _1 , 
which is ~ 50% of their value. We note that our pointing center was offset from that of the 
IRAM observations by (—14", —4"), and Braine et al. (1993) also put a note on the pointing 
uncertainties in the 12 CO(1-0) and 12 CO(2-l) data. Since the 12 CO(1-0) emission is much 
stronger in the SW disk, a strong CO line would show up if the telescope was miss pointed 
toward the SW disk, say at offset (-14," -7"). Comparing with our 12 CO(3-2) and 12 CO(1-0) 
profiles, the IRAM 12 CO(1-0) profile published by Braine et al. (1993) is consistent with the 
profiles in the SW of the nucleus, with the peak shifted to the red side and a broad wing on 
the blue side. Thus it is possible that the IRAM measurement was actually corresponding to 
the SW disk rather than the central nuclear region. This may be the reason for the unusually 
low r 2 i ~ 0.43 reported by BC92. 

Our JCMT 12 CO(2-l) measurements (Figure 3b) yield I 21 = 9.7 ± 1.8 K km s" 1 for 
NGC 157 at offset (0",0"), with a resolution of 21". Convolving the 12 CO(1-0) data to 21" 
resolution we get iio — 9.7 — 9.8 K km s" 1 and therefore r 21 =1.0 in the central 21" region, 
which is consistent with the high excitation suggested by the r 31 ratio towards the same 
position. 
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3.2. The physical conditions of the molecular gas: a study of extremes 

We employ a large velocity gradient (LVG) model (e.g., Goldreich & Kwan 1974) to 
constrain the physical parameters and eventually the CO-to-H2 conversion factor X in the 
extreme ISM environments marked by the two galaxies. This model assumes that the sys- 
tematic motions rather than the local thermal velocities dominate the observed linewidths 
of the molecular clouds. We use this model to fit the observed line ratios for different com- 
binations of (T kin , n(H 2 ), A), where T kin is the kinetic temperature and A = Z co /(dV/dr), 
with Z C o= [ 12 CO /H 2 ] being the fractional abundance of 12 CO with respect to H 2 and dV/dr 
the velocity gradient. The optimum set of parameters is determined by minimizing x 2 - 

Table 4 lists the input r 2 i and r 3i ratios used for our LVG analysis, and the details of 
our method are outlined in Zhu et al. (2003; 2007). The X factor can be directly derived 
from the LVG parameters using the formula 

X = |2^, (3) 

^CO^rad 



where T rad is the radiation temperature for the 12 CO (1-0) line transition. 

In order to reduce the well-known degeneracies in the deduced properties of molecular 
gas from radiative transfer modeling of a few line ratios of the (usually) optically thick 12 CO 
lines we use the constraint of T kin > T dus t (where T dust is derived from FIR and submm data 
in §3.3). This inequality is driven by FUV-heating of gas and dust (usually not thermally 
coupled unless n(H 2 ) > 10 4 cm -3 ), the turbulent heating of the gas (but not the dust), and 
the much more efficient dust cooling (via continuum emission) with respect to that of the gas 
(through spectral line emission). As shown in §3.3, T dust is 32 K for NGC3310 and 23-28 K 
for NGC 157 nuclear region. Hence we estimate T^m = 35-50 K for NGC3310, and 25-40 K 
for the central region of NGC 157. 

A further constraint can be set by adopting only the LVG solutions that correspond to 
gas motions that yield a restricted range of K vw values where 



K ■ 

vir 



(dV/dr] 



obs 



(dV/dr), 



1.54 



Z 



CO 



y/aA 



n(H 2 



10 3 cm 



-1/2 



(4) 



must be K v[r > 1, with K vir = 1 corresponding to the virialized motions of fully self- 
gravitating gas "cells" (Papadopoulos & Seaquist 1999) (a ~ 1 — 2.5, depending on the 
assumed cloud density profile). From the modeling of molecular line emission from quiescent 
to star-forming environments (and using several 12 CO and 13 CO lines) it was found that 
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K vir ~ 1 — 15 (e.g. Papadopoulos & Seaquist 1999), which we adopt here as the permitted 
range of values. 

3.2.1. NGC3310 

The high r 2 i and r 31 ratio in this galaxy can only be fitted with a gas phase that 
corresponds to an optically thin 12 CO(1-0) line. Tellingly no fits are possible for T kin < 30 
K and n < 3000 cm -3 under the constraint K vir < 15. Similar constraints have also been 
reported by Yao et al. (2003) when modeling equally high r 3 i ratios found for galaxies in 
their sample. 

The best fits are found in the range of T kin = 40 - 60 K, n(H 2 )= 5000 - 8000 cm" 3 , 
A = 1.0— 2. 7x 10~ 6 pc [kms -1 ] -1 , which can fit the r 21 and r 31 ratios within the observational 
uncertainties (x 2 < 2). For this range the X factor is 1.3 — 2.6 x 10 19 /(^^) cm -2 [K 
kms" 1 ] -1 . If we further require K vir to be minimum, i.e. a velocity gradient not much higher 
than its virial value, the best fit would then be n=7800 cm -3 , Tki n =50 K and A = 2.7 x 10~ 6 
pc [kins" 1 ]- 1 , with X = 2.6 x 10 19 /(|^) cm" 2 [K kms- 1 ]- 1 ^^ 5). 

A major uncertainty is associated with the adopted CO abundance Zco, usually in the 
range 0.5 — 2.7 x 10~ 4 (Blake et al. 1987; Hartquist et al. 1998) in star-forming clouds of our 
Galaxy. In the case of NGC3310 Pastoriza et al. (1993) found metal abundances 2.5 to 5 
times lower than Solar in its circum-nuclear H II regions, and thus the CO abundance could 
be a factor of 5 lower than the average value of our Galaxy (10~ 4 ). For Z C o = 2 x 10" 5 , the 
X factor would be 1.3 x 10 20 cm -2 [K km s -1 ] -1 (Table 5), which is a factor of 1.5 times lower 
than the Galactic value X G = 2.0 x 10 20 cm~ 2 [K km s^ 1 ]" 1 (e.g. Strong & Mattox 1996), 
with the excitation bias partial offset by the metallicity effect. 

In all the above LVG solutions, the opacity at 12 CO(1-0) is r 10 < 1, quite unlike in typical 
Giant Molecular Clouds (GMC) where this line is optically thick (r 10 > 5). Interestingly, 
only molecular gas in the vicinity of an SNR (and NGC3310 is expected to have lots of 
them) presents high CO ratios (Bolatto et al. 2003) and optically thin CO emission. Even 
in extreme starbursts such as (Ultra) Luminous Infrared Galaxies (LIRGs and ULIRGs) as 
Arp 220 and NGC 6240, r 21 and r 31 ratios are all < 1 (Greve et al. 2009). Thus the molecular 
gas in NGC 3310 seems to be indeed defining the uppermost envelope of average excitation 
conditions known for IR luminous galaxies in the local Universe. 
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3.2.2. NGC157 

The ratios of integrated line intensities for NGC 157 span a wide range, but still within 
that seen in other galaxies (e.g. Yao et al. 2003), making their LVG fitting more straight- 
forward. At the nuclear region conditions appear to be similar to those of NGC 3310 with 
T kin = 30-40 K, n(H 2 )= 2000-5000 cm" 3 , and A = 2.7-7 x lO" 6 pc [km s" 1 ]" 1 (Table 5). 
The range of kinetic temperature is that of the dust since, with limited heating sources in 
the central region of this galaxy (e.g. not much star formation or AGN activity), we don't 
expect the gas temperature to be much higher than that of dust. 

Ratios r 31 < 0.4 such as that found at the offset (—7", —42") where r 31 =0.19, indicates 
that the gas is either very cold (~ 10 K), or of low density (less than 1000 cm" 3 ). For 
r 31 =0.2-0.4 typical LVG fits yield T kin =10-20 K, n = 1000 - 2500 cm" 3 , A = 6 - 12 x 10~ 6 
pc [kms -1 ]" 1 , and X = 2 — 3 x 10 20 cm -2 [K kms -1 ] -1 , while K vir is in the range of 3-4, 
i.e. close to the virial value and thus consistent with what is expected for typical ensembles 
of cold, self-gravitating molecular clouds found in the Galaxy. In general, without internal 
heating sources, (e.g. no star formation) typical gas temperatures in molecular clouds are 
~ 10 — 15 K, set by the equilibrium reached between heating by ambient interstellar radiation 
(G ~ 1 — 5 in Habbing units) and comic rays and C + cooling, while they usually have 
Kvir ~ 1 (to within a factor of a few). 

On the other hand, if we were to set T^ n = 30 K, the best fit would yield n(H 2 ) ~ 1000 
cm" 3 , and A = 2 - 9 x 10~ 6 pc[kms^]^, which gives X = 2.4 - 4.2 x 10 19 cm" 2 [K 
kms" 1 ] -1 . In this case, K vir is in the range 9-30, too high for a quiescent galaxy similar to 
the Milky Way. In §3.3.3, we will see that variations in the SCUBA S4 50 /S 850 flux density 
ratios indicate a variation in the dust temperature, with T dust much higher in the central 
region than in the outer disk. Moreover the global spectral energy distribution (SED), with 
more data points available, is best fitted by a two-component dust model with Tw = 27.9 K 
and Tq = 10.6 K for the warm and the cold dust component respectively. The warm dust 
component is then expected to be associated with the molecular gas in the nuclear region, 
while the cold dust with T c = 10 — 15 K with the gas in the disk and spiral arms. 

3.3. NGC 3310 and NGC 157: The dust continuum emission 

3.3.1. Dust temperature and mass: the methods 

Using the data obtained from the JCMT archive, Stevens et al. (2005) published the 
SCUBA maps of 14 galaxies, including NGC 3310 and NGC 157. They fitted their SEDs with 
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a two-temperature dust model. Their results show that NGC 3310 has the largest fraction 
of cold dust among the 14 sample galaxies. However, there appeared to be a mistake in 
the SED model of NGC 3310 presented in Figure 17 of Stevens et al. (2005), in which the 
submm fluxes were set to an order of a magnitude larger than the real ones, and thus their 
reported dust masses and temperatures are in error. 

We re-analyzed the SCUBA maps, using our data obtained in 2000 (§2.2) and also 
included archival JCMT observations of these galaxies (1999 August 08, 2000 January 19, 
2000 April 13 and 2003 May 20; weather grade 2, 3, 1 and 1 respectively). The total fluxes 
derived (Table 6) are in good agreement with those of Stevens et al. (2005). We then correct 
for 12 CO(3-2) line contamination of the 850 /im fluxes using the formula Sqo = 0.53Jco 
mJy beam -1 (K km s -1 )" 1 (Zhu et al. 2007). This yields a ~ 15% flux correction for 
NGC 3310 (central 30") and a 8-20% flux correction for NGC 157. We furthermore corrected 
the fluxes for contamination of, predominantly thermal, radio emission. Radio continuum 
observations at 5 different frequencies between 57 MHz and 4.85 GHz were available from 
NED for NGC 3310 and data points at 1.4 and 5 GHz for NGC 157. For NGC 3310 we fitted a 
combination of a synchrotron and thermal radio spectrum to the observed radio data points 
by varying the synchrotron spectral index and the synchrotron and thermal flux densities at 
1 GHz. The thermal radio emission has a fixed spectral index of —0.1. Fits were possible 
with a synchrotron spectral index between —0.7 and —0.8 and yielded extrapolated radio 
flux densities between 29 and 48 mJy at 850 /im and between 25 and 44 mJy at 450 /im. 
We adopted the mean value as a realistic estimate. For NGC 157 a detailed fit of the radio 
spectrum was not possible. We estimated the range of the radio continuum contaminations 
by fitting a combined synchrotron and thermal radio spectrum with the steepest (—1.1) and 
the flattest (—0.6) synchrotron spectral index compatible with the radio data points at 1.4 
and 5 GHz and found the range to lie between 5 and 28 mJy at 850 /im and between 3.5 and 
26 mJy at 450 /im. As a realistic estimate for the radio continuum contamination we used a 
synchrotron spectral index of —0.9 and derived radio flux densities of 18 and 16 mJy at 850 
/im and 450 /im. For both galaxies we found that the conclusions of our study are not affected 
by the uncertainty in the radio continuum contamination of our 850 /im and 450 /im flux 
densities. The observed and corrected fluxes (for 12 CO(3-2) and radio contamination) are 
listed in Table 6. 

We then fitted the dust SEDs with two models: (i) a classic two-temperature dust model 
(hereafter 2T model) similar to that used in Stevens et al. (2005) and (ii) a dust grain model 
of Desert et al. (1990, hereafter DBP90 model) which contains Large Grains (LG, (3=2), Very 
Small Grains (VSG, (3=1) and Polycyclic Aromatic Hydrocarbons (PAH). The 2T model uses 
two grey-bodies to account for the warm and cold dust components of temperatures T\y and 
Tc respectively and, in our case, was applied to the data with wavelengths > 60 /im. Both 
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components of the 2T model have the same dust emissivity (3 = 2 in \ v oc (Dunne & 
Eales 2001). The DBP90 model was previously applied to NGC 1569 (Lisenfeld et al. 2002) 
and here we use the same fitting procedure. We adopt the values for the size distribution 
of the different grains for which DBP90 achieved the best-fit for the solar neighbourhood, 
except for the relative abundance of the various grain components which we leave as a free 
parameter. 

The use of the DBP90 model was motivated by the inability of the typical 2-component 
grey-body fit to reproduce the submm "excess" observed in NGC 3310 (see §3.3.2) without 
yielding an unreasonably cold dust temperature and a very large associated dust mass. In 
Appendix A we show how we derive the interstellar radiation field (ISRF) for both galaxies. 
The ISRF in NGC 3310 is much stronger than in NGC 157, even higher than for the star- 
bursting dwarf galaxy NGC 1569 where indications for a boosted VSG population have been 
found (Lisenfeld et al. 2002, Galliano et al. 2003). 

The dust mass M dust in the 2T model was calculated as the sum of the warm (M Wjdust ) 
and the cold (M c dust ) dust using 

S'x , D 2 

w / c ^ 

w/c ^ v ( -^w/ c,dust ) 

where A w / C is the reference wavelength for each component taken as 100 fim for the warm dust 
and 850 /zm for the cold dust and T w / C) d us t the temperature of the warm and cold component, 
respectively. B U (T W / Cj d U st) is the Planck function at 100 and 850 /im, respectively, S' x the 
flux density at 100 and 850 /zm (after correcting the flux at 850 /iin for 12 CO(3-2) and radio 
continuum contamination) and ka w/c the dust opacity per unit dust mass, assumed to be 
0.77 cmV 1 at 850 pim (c.f., James et al. 2002; Dunne & Eales 2001) and 25 cmV 1 at 100 
/iin (Hildebrand 1983). In DBP90, the dust temperature of the LGs was derived by fitting 
equation (5) to the model. We then derived the mass of LG component from the flux at 
850 ^m using the same equation. The masses of the VSGs and PAHs were derived from the 
input mass extinction coefficients, obtained by scaling the values in Table 2 of DBP90. The 
results are presented in Table 7 and Table 8 for the 2T and the DBP90 model respectively. 



M w/ , dust (M ) 



3.3.2. The dust properties, a study of extremes 

The dust SED of NGC 3310 is very different compared to a typical dust SED of a spiral 
(see e.g. Dunne & Eales 2001). The unusual aspect of it is an apparent "excess" emission 
at 850 /im (see Figure 4a), which manifests itself as a relatively shallow submm slope (~ 
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A" 2 - 5 ), derived from the S^o/Ssso ratio of 4.9 (after correcting for CO and thermal radio 
contribution). This is similar to that of the starburst dwarf galaxy NGC1569, which has a 
<S , 45o/<S , 85o ratio of 4.4 (Lisenfeld et al. 2002). By comparison, NGC 157 has 5 450 / 5*850 =6.0 
(after CO and thermal radio correction), and all the galaxies in the SCUBA Local Universe 
Galaxy Survey (SLUGS) have 5 4 5o/5s5o> 6 (Dunne et al. 2000). The value of £450/5850 — 4.9 
for NGC 3310 corresponds to (5 < 2, and constitutes strong evidence of large-scale differences 
of dust properties in starburst environments. This has been suggested earlier using 1.3 mm 
imaging (Dumke, Krause, & Wielebinski 2004), and along with the very high molecular 
gas excitation implied by the high CO line ratios, points to different ISM properties over 
large scales in starbursts. The latter, apart from being interesting in its own right, can 
have serious implications for interpreting molecular line and dust continuum data from the 
numerous starbursts that have now been identified in the distant Universe. 

The need for a different dust SED model in NGC 3310 is exemplified by the unlikely 
results yielded by the classic 2T model (Figure 4a) in which a temperature of 5.7 K is derived 
for the cold component which would then contain the bulk of the mass with Md us t = 2.2 x 10 8 
M (see Table 7), resulting in an unrealistically low GD mass ratio of about 5 within the 
inner 60", the region where most of the dust emission originates (see Table 8 for the gas 
mass). The highest cold dust temperature which is compatible with the error bars of our 
data is 11.8 K. In this case the dust mass would be M dust = 2.8 x 10 7 M and the GD ratio 
= 43, still a low value. The GD ratio for this likely metal-poor object ought to be higher 
not lower than the Milky- Way value of about 100-150. In principle, an underestimate of 
the gas mass could be responsible for the low GD ratio. This has been suggested for the 
Large Magellanic Clouds (Bernard et al. 2008) where cold HI and/or a low density H2 phase 
without CO might cause a low GD ratio. We cannot exclude this possibility in NGC 3310, 
but there are further, stronger arguments against the presence of very cold dust. A first 
argument is that the extremely low dust temperature is also not consistent with the results 
of our CO line excitation analysis. Based on the multi-transition CO data, NGC 3310 has 
a higher gas excitation, and a higher kinetic temperature (T\, in > 30 K) in the molecular 
gas content than NGC 157 (T kin = 10 — 20 K). Thus, we would also expect a higher average 
dust temperature, of the order 30 K or higher, in NGC 3310 than in NGC 157, which is an 
expected difference in the ISM properties of an intense starbursts such as NGC3310. 

A further argument against a very cold dust component is posed by the question where 
the cold dust could hide. Even though such low dust temperatures can be achieved in 
very shielded environments, the dust capable of producing this shielding would emit at 
much higher temperatures. Lisenfeld et al. (2002) estimated that in the case NGC 1569 a 
considerable amount of (warm) dust is necessary to shield the cold dust, so that an extreme 
dominance in mass of the cold dust is impossible to happen. Fischera & Dopita (2008) 
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modeled the SED of dust emitted from self-gravitating, spherical, interstellar clouds which 
allows to base such an estimate on a quantitative foundation. The dust properties that they 
assumed are consistent with Galactic dust. They showed that massive clouds, close to the 
collapsing limit, in a high pressure environment reach the highest opacities. Although in such 
clouds the temperatures of the grains in the inner part can be below 10 K (see their Figure 
8), the emission of this cold dust component does not lead to a flat submillimeter spectrum 
(see their Figures 10 and 11) since the overall emission of the cloud is still dominated by the 
significant mass of warm grains from the outer regions (needed to shield the inner ones so 
that their temperature can drop that low). Thus, very cold dust, even when present in the 
centers of some clouds, is not expected to dominate the dust emission for an ISRF similar 
or stronger than the Galactic one. 

Using the DBP90 model for NGC 3310 on the other hand gives much more reasonable 
results (see Figure 5a), with a LG dust temperature of 32 K, but with an enhanced VSG 
abundance (factor 3.2) and a slightly lower PAH abundance (by a factor of 1.4) compared to 
the Galaxy (see Table 8). The enhancement of VSGs is similar to that found for the dwarf 
starburst galaxy NGC 1569 (Lisenfeld et al. 2002) which is not unexpected given the fact 
that both galaxies are low-metallicity starbursts. A submm excess and enhancement of VSG 
have also been found for other starbursting dwarf galaxies (Galliano et al. 2005). The high 
abundance of VSGs may be the result of grain fragmentation due to numerous SN-induced 
shocks taking place in the ISM. The total dust mass estimated with this model is 2.7 x 10 6 
M Q , with the largest contribution coming from the LG component (2.1 x 10 6 M ). The 
average dust temperature deduced for the LGs is 32 K and is consistent with that of the H 2 
gas. The GD mass ratio within the central 60" (the region where most of the dust emission 
emanates from) is 461. This value is ~ 3 times higher than the Galactic value, but within 
the range expected for the low metallicity reported for NGC 3310 (see §3.4.1). 

The good fit to the shallow submm SED of NGC 3310 is partly due to the fact that 
the stochastically heated VSGs possess a broad temperature distribution going down to very 
low values. Most importantly, however, the low dependence of the extinction coefficient on 
the wavelength ((3 = 1) of the VSG component in the DBP90 model is responsible for the 
good fit. The empirical nature of the DBP90 leaves some uncertainty about whether this 
low value of (3 is correct. In any case, a firm conclusion from the good fit of the DBP90 
model together with the impossibility to use a two-temperature model (with (5 — 2) is that 
we need a low value of (3 in order to fit the shallow submm spectrum of NGC 3310. 

Low values of (5 are expected for amorphous carbon (Koike et al 1980), aggregates of 
silicates and graphite in a porous structure (Mathis & Whiffen 1989) as well as for small amor- 
phous grains (Seki & Yamamoto 1989, Tielens & Allamandola 1987). In general, amorphous 
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rather than crystalline structure and increase of the grain size (e.g. Mannings & Emer- 
son 1994) are the main reasons for expecting /3 < 2. Furthermore, laboratory experiments 
have shown that f3 can be temperature dependent and decreases at very low temperatures 
(Agladze et al. 1996) which could indicate that dust in the inner parts of dense clouds might 
show different properties than dust outside. Going into a similar direction, Meny et. al. 
(2007) proposed a dust model based on physical properties of amorphous solids which ex- 
hibits a broad emission spectrum whose detailed properties dependent on temperature and 
which might explain the submm excess. Thus, an alternative explanation for the flat submm 
SED could be the change of dust properties with temperature. 

The SED of NGC 157 is similar to that of our Galaxy and both dust SED models can 
fit the data well, but with a dust mass different by a factor of 3-4. From the 2T model 
we derived dust temperatures of 28 K and 11 K for the warm and the cold dust component 
respectively (see Figure 4b and Table 7), in good agreement with the results of Stevens et 
al. (2005). The dust mass is 2.34 x 10 8 M , corresponding to a global GD mass ratio of 124, 
similar to the Galactic value. From the DBP90 model (Figure 5b) the temperature derived 
for the large grains is 22.5 K, with a total mass of 6.2 x 10 7 M Q (Table 8), and we achieved a 
good fit for VSGs and PAH abundances similar to those in the Galaxy. This model yields a 
GD mass ratio of 468, which is 3 times higher than the Galactic value. This is mostly due 
to fitting the LGs with a single temperature component of 22.5 K. As shown in §3.3.3, the 
flux density ratios of S450/S850 suggest that the dust temperature indeed varies by a factor 
of 2 across the disk of NGC 157, thus the two-temperature model may better reflect the dust 
distribution in normal spiral galaxies such as NGC 157. 



3.3.3. The S^so/Ssso ratio: dust temperature variations, and VSGs versus LGs 

Having determined the proper dust model for each galaxy we can now examine whether 
there are any discernible variations of the dust temperature within these two systems. Such 
variations have been identified via submm imaging in the starburst/AGN NGC 1068 by 
Papadopoulos & Seaquist (1999) using the S i50 /Ss5o ratio. This ratio can be a good indicator 
for temperatures of Td us t < 30K, which characterize quiescent H2 or H I-rich regions with 
little ongoing star formation. For a single dominant dust component this ratio can be 
expressed as 

= l-^ Wl )» (6) 
with (3 being the emissivity index (2 in our case). For a quiescent spiral such as NGC 157, 
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the S'450/5'850 variation reflects a change in the dust temperature. In Figure 6 we show the 
azimuthally averaged S450/S850 ratio (after smoothing the 450/zm map to the resolution of 
the 850/^m map; Figure 6a) and dust temperature profile as a function of radius (Figure 6b) 
for both galaxies. We don't expect that equation (6) can be applied to NGC 3310, given that 
in this galaxy: a) the large grains have Td us t > 30 K (and thus their submm emission lies in 
the Rayleigh Jeans domain), and b) the VSG component could be present throughout the 
compact submm emitting region of this galaxy. To demonstrate this, Figure 6b shows that 
the "dust temperature" of NGC 3310 derived with equation (6) is unreasonably low, much 
lower than that of NGC 157. Therefore, the variation of the S450/S850 ratio for NGC 3310 
in Figure 6a, may not reflect a simple variation in T dust , but variations of the "effective" 
submm-deduced j3, possibly due to a change of the VSG/LG mass ratio, something that our 
global dust SED fits of this system cannot discern. We will address this issue in a future 
paper, with detailed analysis of MIPS data from Spitzer. 

In NGC 157, the dust temperature changes from 23-28 K in the central region to less 
than 15 K in the outer disk (Figure 6b). This is consistent with the variation in the CO 
gas excitation temperature estimated from the r 31 ratio, and again suggests that a two- 
temperature dust model represents a good fit to the global SED in this galaxy. 

3.4. CO, H I, and dust distributions: an empirical test of the X factor 

In this section we study the variation of the GD mass ratio in NGC 3310 and NGC 157 
in order to explore whether we can assume a constant GD mass ratio for the H I-dominated 
and CO/H 2 -dominated regions to get an independent, empirical test of the X factor (for a 
similar method applied in the Milky Way see Reach et al. 1998). The dust mass is estimated 
using equation (5) and the H I and H 2 mass can be estimated using the following formulae: 



where D is the distance to the galaxy in Mpc, Shi and Sqo is the H I and 12 CO(1-0) 
integrated flux in Jy km s -1 . Equation (8) for the Mh 2 estimate is essentially equation (4) in 
Braine et al. (2001) when using the Galactic X factor X G = 2.0 x 10 20 cm -2 (K km s -1 ) -1 or 
4.3 M Q pc~ 2 (K km s -1 ) -1 (including a 36% mass correction for helium). Our LVG analysis 
shows that the factor X/X G is approximately unity in the disk of NGC 157, but could be 
0.65 in NGC 3310 and 0.28 in the nuclear region of NGC 157. 



M m (M ) = 2.36 x 10 5 £> 2 Shi 
M H2 (M ) = 1.073 x 10 4 (X/X G )D 2 S CO , 



(7) 
(8) 
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For both galaxies, the 12 CO(3-2) maps are more extended than those of 12 CO(1-0) and 
cover the majority of the CO distribution. They have higher S/N ratios and hence we used 
them to represent the CO distribution, bearing in mind the variation of r 31 . Also, we used 
the uncorrected SCUBA image to get the highest S/N in order to trace the dust emission to 
the largest extent. The CO contribution is on average about 15% and thus not significant. 

34.1. NGC3310 

Figure 7a shows the 12 CO(3-2) contours overlaid on the SCUBA 850 /im image while 
Figures 7b and 7c show the SCUBA 850 /im and 12 CO(3-2) contours respectively overlaid 
on the H I map at 20"resolution, taken with the WSRT interferometer (Kregel & Sancisi 
2001). The 12 CO(3-2) map covers the majority of the CO emission in the central 50" x 50" 
region, while the dust emission can be detected out to a 60" x 90" region, thanks to the high 
sensitivity of SCUBA. The IRAM 12 CO(1-0) and 12 CO(2-l) maps published by Mulder et 
al. (1995) indicated the existence of weak CO emission at ±40" off the center, slightly larger 
than the area we detected. The total 12 CO(1-0) flux from the 45-m map is 101 Jykms -1 . 
Comparing to the total flux of 140 Jykms -1 reported by Young et al. (1989), our map 
covers 72% of the CO, mostly in the central region. 

The H I emission is much more extended than the CO and SCUBA 850 /zm emission 
and it dominates the gas content outside 60". Unlike NGC 157 (see §3.4.2), there are no 
regions that are dominated by CO. In the central 20", the Mhi/Mh 2 ratio is about 0.8 if 
X = 1.3 x 10 20 cm~ 2 [K kms _1 ] _1 is assumed. However, the H I map at 20" resolution is 
missing fluxes by about 30%, which was estimated by comparing to the 60"map (also from 
Kregel & Sancisi 2001), after convolving to the same resolution. Thus we apply a factor of 
1.3 to the high resolution H I map, and the resulting Mm/ Mh 2 is revised to 1.03, i.e. neither 
H I nor H 2 dominates the gas content in the nuclear region. The 850 /im flux is 67 mJy 
in the central 20", and we estimate the dust mass by scaling this flux with the total flux 
to obtain the fraction of M dust , derived from the DBP90 model, in this region. The final 
derived GD mass ratio is 312 (Table 9). If a Galactic X factor is used, the central region 
would have 54% more molecular gas. The M H i/M H2 and GD mass ratio would become 0.67 
and 395, respectively. 

Similarly, in the central 60" region (r < 30"), the GD mass ratio would be 461 (the 
detailed numbers are listed in Table 9). Our CO data does not cover the entire 60" region, 
but Mh 2 contributes only about 27% to the total gas mass, thus a 20-30% under-estimate 
in the CO fluxes would have a limited effect on the GD mass ratio. 
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In the outer region with r > 40" where little CO emission is found, there is some weak 
850/im emission associated with the H I. To improve the S/N ratio, we averaged the S$5o 
fluxes azimuthally in the outer disk and used the azimuthally averaged Shi/iS^o ratio to 
estimate the GD mass ratio which is essentially the M H i/M dust ratio as: 

M m /M dnst = 593(e 16 - 8 / Td - - l)" 1 . 

Using the value T dust = 32 K from the DBP90 model, we estimate a GD mass ratio of 859. 
If the VSG component is not enhanced in the outer disk (possibly due to less star formation 
activities) and T dust drops to 20 K, the GD mass ratio would be 456. 

In all cases, the resulting GD mass ratio in NGC3310 is much higher than the Galactic 
value. The high GD mass ratio could be related to the low metallicity in the galaxy. Pastoriza 
et al. (1993) found metal abundances 2.5 to 5 times lower than Solar in circumnuclear H II 
regions of NGC 3310. It is thus possible that this galaxy has less dust compared to a solar- 
metallicity galaxy. The lack of dust is consistent with the low global M H2 /M H i ratio in 
NGC 3310 as molecules are formed on the surfaces of dust grains. 

The GD mass ratio in NGC 3310 shows larger variations than in NGC 157 (see §3.4.2). 
Also, due to the uncertainties in estimating the fraction of the VSG component and hence 
the dust temperature, the derived local GD mass ratio has a large uncertainty, by a factor 
of at least 2. Combined with the fact that there are no H 2 dominated regions in this galaxy, 
we conclude that it is not practical to use the GD mass ratio to derive the X factor in 
NGC 3310. 



34.2. NGC157 

Figure 8a shows the 12 CO(3-2) integrated intensity contours overlaid on a SCUBA 850 
/im image. The CO distribution is in general consistent with that of the dust emission. A 
major difference is visible in the nucleus of the galaxy, where 12 CO(3-2) emission is strongest 
but the dust emission is relatively weak. In fact, the strongest 850 fim emission is located 
in the southern spiral arm, close to the 12 CO(1-0) peak (see Figure 2). This difference is 
probably due to an excitation effect. The 12 CO (3-2) is tracing hot molecular gas and its 
distribution is more consistent with that of 450 fim emission (both have a strong emission 
near the centre). This suggests that different dust temperatures and X factors (for different 
excitations) should be used to derive the M H2 /M dust ratios in different regions. 

Figures 8b and 8c show the SCUBA 850 //m and 12 CO(3-2) contours respectively over- 
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laid on the H I map at a resolution of 18" x 12", obtained by Ryder et al. (1998) using the 
VLA. The H I map was constructed using the uniform weighting schemes (referred as UN 
image) which has a high resolution and shows the structure of the inner region, but it misses 
some extended structure hence is not good for deriving the H I column density. We use it to 
compare the H I morphology with other ISM components in the inner region. Ryder et al. 
(1998) also published a low resolution (41.5" x 27.2") VLA map processed with the nature 
weighting schemes (referred as NA image), which is expected to recover all the H I fluxes. 
We have converted the UN image to 42" x 27" resolution to compare with the NA image and 
found that about 40-50% of fluxes could be missed in the UN image. 

Based on these maps we divide the galaxy into 3 parts: (i) the H 2 dominated region in 
the central 21", (ii) the H I dominated outer disk with r > 50" (iii) The regions between (i) 
and (ii). In the following we discuss them separately. 

(i) In the central 21" region, the Mhi/Mh 2 ratio is only 0.2 based on the H I fluxes from 
the high resolution UN image, after applying a correction factor of 2 for missing fluxes. The 
central region is clearly dominated by molecular gas and thus the uncertainty in the H I 
mass estimate has a limited effect on the total GD mass ratio. Line 4 of Table 9 lists the 
mass ratios for the central 21". For dust mass estimates we used T dU st=23 K which is the 
temperature we derived for the central region of this galaxy (§3.3.3). The GD mass ratio 
turns out to be 54 if a low X factor (5.6 x 10 19 cm -2 [K kms -1 ] -1 ) is used. 

(ii) In the outer disk with r > 50", most of the gas are in the atomic phase, and the 
GD ratio is effectively M H i/M dust . Similar to the case of NGC3310 outer disk, we used the 
azimuthally averaged Shi/Ssso ratio to derive the Mm/Mdust ratio as: 

M ffl /M dust = 95(e 16 - 8 / Td - - l)- 1 

If Td U st = 15 K, as derived from equation (6) in §3.3.3, the GD mass ratio would be 46. This 
is consistent with the GD mass ratio in the central H 2 dominated region. 

(iii) For the region with 10" < r < 50", we used the NA image to derive the H I fluxes 
and convolved the CO and 850 /im data to a resolution of 42" x 27". At such a low resolution, 
we select three regions to derive the GD mass ratios: the center, the NE and the SW disk 
that are 40" offset from the nucleus. Table 9 lists the CO, H I and dust mass in these 
three regions and the GD mass ratios. The X factor was assumed to be 1.1 x 10 20 cm -2 [K 
kms -1 ] -1 ^ the center, but 2.0 x 10 20 in the SW and NE spiral arms. The dust temperature 
is 23 K in the center but 15 K at 40" away from the nucleus, according to the temperature 
distribution in Figure 6b. 
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Our result indicates that the variation of the GD mass ratio is less than 40%. Thus 
we can assume a constant GD mass ratio of 62 ± 20 to estimate the X factor in the H 2 
dominated region, which yields X = 6.6 x 10 19 cm' 2 [K kms _1 ] _1 for the central 21", in good 
agreement, and independently, with the results from our LVG analysis of the CO lines. 

4. Discussion 

4.1. NGC 3310 and NGC 157: two extremes in ISM excitation 

Few molecular line ratio surveys of large samples of galaxies exist, and those available 
are of the r 2 i CO line ratio (Braine & Combes 1992; Papadopoulos & Seaquist 1998), while 
smaller samples have been observed in the more excitation sensitive r 31 ratio (Yao et al. 
2003; Narayanan et al. 2005). With the new array receiver HARP and the ongoing JCMT 
nearby galaxies legacy survey (e.g. Wilson et al. 2009), more r 3 i studies of nearby galaxies 
are expected. 

In the Braine & Combes (1992) survey the starburst NGC 3310 and the quiescent spiral 
galaxy NGC 157 have respectively the highest and lowest r 2 \ values ever observed in galaxies. 
Although our observations did not confirmed these extreme values, the r 3 i ratio of 1.5 in 
NGC 3310 remains among the highest of such ratios measured in galaxies, while in NGC 157, 
an extremely low ratio of r 31 = 0.2-0.4 is found in the outer disk. These two galaxies indeed 
represent two extremes in ISM excitation in external galaxies, allowing us to study the range 
of the X factor in extreme environments. 

Our results showed an X factor varying by a factor of ~ 5 between the two galaxies 
and within a galaxy where strong excitation gradients are present. In the disk region of 
NGC 157 the X factor is close to the standard Galactic value found for Galactic GMCs at a 
galactocentric radius of ~ 5 — 10 kpc, while at its center it could be lower by a factor of ~ 5, 
a trend also verified for the Milky Way (Sodroski et al. 1995). This good correspondence 
demonstrates the utility of excitation analysis of multi-transition CO data to yield good con- 
straints on the X factor in external galaxies, even when a significant degree of degeneracy 
in the deduced average physical conditions of the molecular gas remains. The same method 
has been used by Zhu et al. to study the X factor in interacting galaxies, e.g. the Antennae 
(Zhu et al. 2003) and the Taffy galaxies UGC 12914/15 (Zhu et al. 2007). These galaxies 
represent other types of extreme environments for molecular clouds, with strong tidal inter- 
action, direct cloud-cloud collision, and diffuse clouds outside a galactic disk environment 
(in the case of Taffy galaxies). The X factor was found to be in the range of 0.1 - 0.2 Xq 
in the these galaxies. Yao et al. (2003) also found on average X = O.IXq based on their 
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analysis of the r 3 i ratios in the central regions of SLUGS galaxies. 

In most cases studied so far, a lower than standard X factor is associated with non- 
virialized clouds rather than non-standard molecular gas physical conditions. According to 
equations (3) and (4) 

T Tad dV/dr T rad vir! 1 ) 

with K vlT = 1 for ensembles of virialized clouds (the ensemble itself in general does not need 
to be virialized, and its motions ordinarily are dictated by the general galactic potential). 
Most molecular clouds in a quiescent galaxy have K yu ~ 1 and X ~ Xq. Our LVG excitation 
analysis find a vW^~rad varying at most within factors of < 3, while major variations come 
from K V1I > 1. In the Milky Way, very large local velocity gradients (greater than 10 km 
s _1 pc _1 ) with K vir 3> 1 are found near pre-star-forming clouds (Falgarone et al. 1998) as 
well as in diffuse clouds towards Orion (Knapp & Bowers 1988). Falgarone et al. (1998) 
found high velocity gradients associated with macroturbulent gas "cells" of size < 200 AU 
and densities of n(B.2) ~ 10 3 cm -3 , which is consistent with the parameters found in our 
LVG analysis. 

In the nuclear regions of ULIRGs, K vir can be significantly higher than unity, due to a 
large scale "mixing" of molecular gas with significant amounts of stellar mass, resulting in 
X ~ 0.2Xg (Solomon et al. 1997; Downes & Solomon 1998). Actually our method of using 
the LVG radiative transfer models to obtain fT vir values can in such a case (with gas and star 
mixtures) yield the factor f g = M H2 /M tota i = (AV7AU vir ) 2 = K~?, independent of methods 
that model the detailed spatial and velocity distribution of molecular gas in such systems, 
which demand high resolution interferometric imaging (e.g. Downes & Solomon 1998). 

Summarizing our findings for the variations of the X factors in NGC 3310 and NGC 157, 
as well as in the aforementioned ULIRGs and interacting galaxies (Zhu et al. 2003, 2007), 
we find that the X factor is generally in the range of 0.1-0.5 X G in extreme environments, 
which is mostly due to the variations of the parameter K vir ranging from 1 to ~ 15. 



4.2. Molecular gas estimates in extreme environments: cautionary tales 

The discovery of a significant population of dust-enshrouded and CO-bright starbursts 
at high redshifts (e.g. Smail, Ivison, & Blain 1997; Hughes et al. 1998; Greve et al. 2005), 
as well as the discovery of extended outer regions dominated by cold dust (and concomitant 
gas) in several spirals (Nelson, Zaritsky, & Cutri 1998; Xilouris et al. 1999; Papadopoulos 
& Seaquist 1999; Thomas et al. 2001, 2002) give to our study of extreme ISM extragalactic 
environments a wider importance. 
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For high-z galaxies for which multi-transition CO data are rarely available we can only 
make an educated guess of the X factor based on our knowledge of local starbursts, frequently 
used as their low-z counterparts based on the high star formation rates per molecular gas 
mass found for them. A gas dynamics driven by fast-evolving galactic potentials associated 
with mergers and/or frequent SNR shocks could be the decisive similarity among them that 
makes claims about similarly low X factors in both populations credible. Indeed if these low 
X values are driven mostly by large velocity gradients rather than gas excitation differences, 
highly turbulent environments would be expected to be present in them. However, it has 
been shown that interstellar turbulence decays quite rapidly, on timescales of the order of 
the free-fall time of the system (e.g. Klessen et al. 2000; Mac Low et al. 1998). To maintain 
the observed long lifetimes, turbulence in molecular clouds must be constantly driven, either 
externally (e.g. via galaxy interactions) or internally via star formation. Therefore, generally 
a gas-rich galaxy with a high star formation efficiency, indicated by a high Lfir/Lco? may 
have systematically lower X factors. Indeed, the Lpm/^co ratio is much higher in NGC 3310 
and ULIRGs than in quiescent spirals. Tacconi et al. (2008) shows an anti-correlation 
between the X factor and the gas surface density (see their Figure 10), implying that galaxies 
with more dense gas and therefore higher star formation efficiencies would have lower X 
factors. We will explore this effect more quantitatively in a future paper. 

The other method of deriving H 2 gas mass from the dust mass, a benchmarking of the 
GD mass ratio from the H I-rich part of the galaxy, while independent of the X factor and 
thus valuable, cannot be currently applied to high redshifts galaxies. Its use is restricted in 
the local Universe of z < 0.3 where H I can be imaged by current interferometers. More- 
over, the presented evidence for different dust properties in the intense starburst galaxy 
NGC 3310, while extremely interesting by themselves (especially if found generally true for 
similar galaxies) complicates further the task of estimating the H 2 gas mass following this 
method, even if H I, IR/submm dust continuum and CO images are available. 



5. Summary and Conclusions 

We present fully-sampled 12 CO(3-2) and 12 CO(1-0) observations and 850/im and 450/im 
dust continuum imaging, all at a matched beam resolution of HPBW ~ 15", for the starburst 
NGC 3310 and the quiescent spiral NGC 157 that represent two extremes of molecular gas 
and dust excitation conditions in IR-luminous galaxies. Including available sensitive H I 
maps convolved to a similar resolution, we have obtained a detailed view of all phases of 
the neutral ISM in these galaxies, and an examination of the standard method of using 
12 CO(1-0) brightness in deducing H 2 gas mass. Our conclusions are summarized as the 
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following: 

(1) The extremely high r 2 i in NGC3310 and low r 2 \ in NGC 157 reported by Braine et 
al. (1993) are not confirmed by our observations. The r 31 and r 2 i ratios in NGC 3310 are 
still at the high end; while in NGC 157 the r 3 i ratio is found to be low only in the outer disk. 

(2) NGC 157 shows large variations in the r 31 ratio, which is interpreted as variations 
in excitation condition and hence the X factor. The X factor is close to the Galactic value 
X G in the quiescent regions of the outer disk where the temperature is cold (no heating) 
and the velocity gradient is low. In the nuclear region, on the other hand, the temperature 
and gas velocity gradient are high, resulting in an X factor being ~ 4 times lower than Xq, 
assuming Z co is 5 x 1CT 5 . 

(3) The gas-to-dust mass ratio found for NGC 157 is consistent with that of the Galaxy, 
and using its benchmarked value in its H I-rich part (assuming it stays the same over its 
entire disk) independently yields an X factor of ~ 4 times lower than its Galactic value in 
its H 2 -rich nuclear region. 

(4) The X factor in NGC 3310 can be ~ 8 times lower than its Galactic value for typical 
CO abundances Z C o = 10~ 4 . If one assumes that the lower metallicity values deduced for 
a few H II regions in this galaxy are representative of its bulk ISM (Zqo = 2 x 10 -5 ), then 
the X factor is only ~ 2 times smaller than the Galactic values with the excitation effects 
mostly offset by the lower metallicities. 

(5) We deduce an unexpected and substantial contribution from a potentially very dif- 
ferent grain population in the starburst NGC 3310, manifesting itself as an observed submm 
"excess" which a classic 2T model fit to its dust emission SED would otherwise interpret as 
extremely cold dust with Tq ~ 5 — 11 K which corresponds to gas-to-dust mass ratios of 
~ 5 — 43, very unlikely for a starbursting and most likely metal-poor system. We find that 
a low value for the frequency dependence of the extinction coefficient, f3 ~ 1, is required 
to explain the observations. We find that a small grain population, a possible result of the 
extreme radiative and SNR-dominated ISM environment of this extreme system, can account 
for this submm "excess" without postulating very cold dust temperatures in such an unlikely 
environment. We also discuss other possible explanations. 

The authors wish to thank the staff of the JCMT for their generous assistance during an 
extended period of observations. The support personnel at the 45-m telescope in Nobeyama, 
Japan is also gratefully acknowledged for their assistance. We also thank Stuart Ryder and 
Michiel Kregel for providing H I images. UL acknowledges financial support from the research 
project AYA2007-67625-C02-02 from the Spanish Ministerio de Ciencia y Education and 
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from the Junta de Anaducia. In addition, we thank the anonymous referee for stimulating 
suggestions and comments. 

A. The Interstellar Radiation Field of NGC 3310 and NGC 157 

In order to calculate the dust emission, in the concept of the DBP90 model, we need to 
quantify the Interstellar Radiation Field (ISRF) in which the dust is immersed. We derive 
this radiation field from the integrated extinction corrected fluxes in the UV and optical 
wavelength (listed in Table 10), assuming a spherical symmetry of the emitting region with 
a radius of 3.2 kpc for NGC 3310 and 12.7 kpc for NGC 157. Both radii are estimated from 
the UV (GALEX archival data reduced by us), optical and SCUBA maps. The small size 
in NGC 3310 reflects the fact that both the UV and optical radiation and the dust emission 
are very concentrated to the center. Figure 9 shows the ISRF field for both galaxies (see 
the figure caption for explanation of the different symbols). We can see that the ISRF in 
NGC 3310 is even stronger than that of NGC 1569, which could enhance the population of 
very small grains. 
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Fig. 1.— The 12 CO(3-2) profiles (red) overlaid on the 12 CO(1-0) profiles (green) of 
NGC3310. The background grey-scale map is the SCUBA image at 850yum. The spec- 
tra have been smoothed to a spectral resolution of 21.6 kms -1 . The velocity range is 660 to 
1320 kms -1 . In each plot the Y-axis is T m b which ranges from -0.02 to 0.15 K. 
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Fig. 2.— The 12 CO(3-2) profiles (red) overlaid on the 12 CO(1-0) profiles (green) of NGC 157. 
The background grey-scale map is the SCUBA image at 850/zm. The velocity range is 1360 
to 1996 kms _1 for 12 CO(1-0) but 1300 to 2020 kms _1 for 12 CO(3-2). In each plot the Y-axis 
is T m b which ranges from -0.05 to 0.18 K. 
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Fig. 3.— The 12 CO(2-l) spectrum observed with the JCMT (solid line) overlaid on the 
12 CO(1-0) spectrum (dashed line) from the NRO 45-m for NGC3310 (a) and NGC 157 (b). 
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Fig. 4. — Two-temperature dust model fitting to the SEDs of the two galaxies (a) NGC 33 f 
and (b) NGC f 57. Two grey-bodies were fit to the data (for wavelengths > 60 /zm ) to account 
for the warm and the cold dust components. 




Fig. 5. — The SED fitting using the Desert et al. 
NGC 157 (b). 



(1990) dust model for NGC 3310 (a) and 
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Fig. 6. — Azimuthally averaged profiles of (a) the S4 50 /S S50 flux ratio, and (b) the dust 
temperature derived from equation (6) for NGC3310 and NGC 157. A radial step size of 
16", comparable with the 850/im beam size, was used. 
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Fig. 7. — Morphological comparison between the different ISM tracers (H I, 850/im and 
12 CO(3-2)) for NGC3310. In Figure 7a, the 12 CO(3-2) contours are overlaid on top of the 
SCUBA 850 /im map. The contour levels are 4, 5.6, 8, 11.2 and 16 Kkms" 1 . In Figure 7b, 
the 850 /im contours are overlaid on top of the the H I map. The contour levels are 10, 14, 
20 28, 40 and 56 mJy/beam. Figure 7c presents the 12 CO(3-2) contours on top of the H I 
map. The contour levels are as in Figure 7a. In all maps the beam size (15") is indicated on 
the bottom-left corner. 
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Fig. 8. — Morphological comparison between the different ISM tracers (H I, 850/im and 
12 CO(3-2)) for NGC157. In Figure 8a, the 12 CO(3-2) contours are overlaid on top of the 
SCUBA 850 /im map. The contour levels are 4, 5.7, 8, 11.4 and 16 Kkms" 1 . In Figure 8b, 
the 850 /im contours are overlaid on top of the the H I map. The contour levels are 17, 23, 
34 and 46 mJy/beam. Figure 8c presents the 12 CO(3-2) contours on top of the H I map. 
The contour levels are as in Figure 8a. In all maps the beam size (15") is indicated on the 
bottom-left corner. 
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Fig. 9. — The intensity of the interstellar radiation field of NGC 3310 and NGC 157 compared 
to that of the low-metallicity dwarf galaxy NGC 1569 (see Appendix A for details). The open 
symbols (triangle and square) are derived from the observed data listed in Table 10, assuming 
the emitting region as described in Appendix A. The filled symbols give the interpolation to 
the observed data points used as input in the DBP90 model. The crosses give the interpolated 
ISRF of the dwarf starburst galaxy NGC 1569 for comparison. 



Table 1: Calibration data of the sub-mm SCUBA 
observations. 



Wavelength 


#HPBW 


G,SG/G a 


&SKY b 


(H 


(arcsec) 


(%) 


(mJy beam" 1 ) 


850 


15.28 


221, 16 


4.2 


450 


9.28 


332, 32 


29.0 



"The flux gain in units of Jy beam -1 V -1 and its 
fractional uncertainty estimated from an extensive set of 
beam maps on CRL618. 

b The sky a in the final images. 



Table 2: 12 CO(1-0) and 12 CO(3-2) integrated intensities in NGC3310. 
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a Both the /10 and I32 intensities are in main beam temperature. Non-detections 
for 12 CO(1-0) have an upper limit of 3 sigma. 

b The 12 CO (3-2) line intensity has been converted from 14" to 15" to match the 
12 CO(3-2) HPBW. 



Table 3: 12 C0(1 



-0) and 12 CO(3-2) 



integrated 



intensities in NGC 157. 
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a Both the iio and ^32 intensities are in main beam temperature. For non- 
detections, upper limits of 3 sigma are listed. 

b The 12 CO(3-2) line intensity has been converted from 14" to 15" to match the 
12 CO(3-2) HPBW. 



Table 4: Integrated intensity ratios at the central 21 
region. 



Source 


ho 


hi 


I32 


T21 


^31 


T31 




K km s _1 


K km s" 1 


K km s _1 








NGC 157 


9.7 


9.7 


8.34 


1.0 


0.86 


1.28 


NGC3310 


6.62 


9.8 


9.88 


1.48 


1.49 


1.17 



a The r3i ratio at a 21" resolution. 
b The r3i ratio at a 15" resolution. 



Table 5: LVG Analysis Results. 



Source 


n(H 2 ) a 


1 kin 


Z co /{dV/dr) a 


no " 


X c 


K ■ d 


x 2e 




10 3 cm -3 


K 


10" 6 pc (km/s)" 1 




1Q 19 






NGC 3310 


7.8 


50 


2.7 


0.3 


13.0 


10.1 


0.9 




5.0-8.0 


40-60 


1.0-2.7 


0.18-0.75 


6.3-13.0 


10-15 


<2 


NGC 157C 


3.5 


30 


2.7 


0.7 


4.6 


6.1 


0.6 




2.0- 5.0 


30-40 


2.7- 7.0 


1.3-1.6 


4.6-6.1 


6-13 


<2 


NGC 157S 


1.0-2.0 


10-20 


6-12 


1.5-6.7 


20-30 


3-4 


<2 



a The LVG parameters described in §3.2. 

6 The optical depth at the 12 CO(1-0) transition. 

c The X factor in the unit of cm -2 (Kkms -1 ) -1 , assuming Zqq, = 5 x 10 -5 
in NGC 157 and Z CO = 2 x 10~ 5 in NGC 3310. 
d The parameter K vu . defined in equation (4). 
e The value of x 2 associated with the fit. 

Notes. NGC 157C and NGC 157S correspond to the central region and and 
the southern disk of NGC 157 respectively. For each position, the first row is 
the best fit parameters and the second row is the estimated range within the 
observed uncertainty. For NGC 157S the r%\ ratio is in the range 0.2 - 0.4 and 
thus a range of values is given for the parameters. 



Table 6: Far-IR and sub-mm flux densities for NGC3310 and NGC 157. 









Integrated flux densities 












NGC 3310 






NGC 157 




A 


observed 


CO radio 


dust 


observed 


CO 


radio 


dust 


/im 


Jy 


Jy Jy 


Jy 


Jy 


Jy 


Jy 


Jy 


12 


1.54± 0.03 




1.54± 0.15 


1.61± 0.04 






1.61± 0.16 


25 


5.32 ±0.03 




5.32 ±0.53 


2.17± 0.04 






2.17 ± 0.22 


60 


34.6 ±0.05 




34.6 ±3.46 


17.9 ± 0.05 






17.9 ± 1.79 


100 


44.2 ±0.17 




44.2 ±4.42 


43. 4± 0.1 






43.4± 4.34 


450 


1.32 ±0.26 


0.035 


1.29±0.26 


4.68± 0.93 




0.018 


4.66 ±0.93 


850 


0.36 ±0.07 


0.051 0.039 


0.27 ± 0.07 


0.93± 0.19 


0.14 


0.016 


0.77±0.19 







Notes. The flux densities at 12, 25, 60 and 100 /im are from IRAS (Soifer et 
al. 2003). The errors for the observed flux densities represent la uncertainties 
of the data, whereas we adopt a 10% error as a more realistic value for the dust 
flux density at these wavelengths. The fluxes at 450 /im and 850 fim are derived 
by us by the integration over the total extent of the corresponding maps from 
the JCMT archive. In order to derive the flux density corresponding to the 
pure dust emission (column "dust" ) we have to correct for the contamination 
of the 12 CO (3-2) at 850 |im (column "CO") and the radio contamination 
(column "radio") at sub-mm wavelengths (see §3.1.1). The corresponding 
values are listed. 



Table 7: SED fitting result from the 2T model. 



Galaxy 


Tw 


T c 


-^w,dust 


-^c,dust 


M dust 




K 


K 


10 s M 


10 8 Mq 


10 8 M 


NGC 3310 


34.4 


5.8 


0.05 


2.19 


2.24 


NGC 157 


27.9 


10.6 


0.18 


2.34 


2.52 



Table 8: SED fitting result from the DBP90 model. 



Galaxy 


^dust 


M LG 


Mvsg 


Mpah 




K 


10 6 M 


10 6 M 


10 6 M 


NGC3310 


32.0 


2.1 


0.5 


0.1 


NGC157 


22.5 


54 


4.5 


4.1 



Table 9: Mass of ISM components and the gas-to-dust ratios 



Position" 


Shi 


5 'i2co(i-o) 6 


M m 


M(H 2 ) c 


X d 


5*850 b 


M dust e 


J dust 


M H2 +M m 
M dust 




Jy km s -1 


Jy km s _1 


10 s M 


10 8 M 


10 19 


Jy 


10 7 M 


K 


NGC 157 




















Centre 


1.40 


71.3 


4.1 


5.16 


11 


0.091 


1.15 


23 


79 


SW 


2.28 


128. 


6.6 


16.9 


20 


0.125 


3.1 


15 


76 


NE 


2.03 


97.9 


5.9 


12.9 


20 


0.106 


2.6 


15 


72 


21" 


0.11 


45.6 


0.33 


1.68 


5.6 


0.029 


0.37 


23 


54 


Total 


78 


500 


225 


65.7 


20 


0.8 


23.4 


10.6 


124 


NGC 3310 




















20" 


0.83 


28.2 


0.92 


0.89 


13 


0.067 


0.058 


VSG 


312 


60" 


8.2 


101 


8.8 


3.2 


13 


0.29 


0.26 


VSG 


461 


Total 


69 


140 


73.9 


4.4 


13 


0.31 


> 0.27 f 


VSG 


<2900 



a The first three positions of NGC 157 are in a 42" x 27" region (resolution of the 
VLA HI data with nature weighting) at the center and two spiral arms 40" offset 
from the center. Line 4 is for the central 21" region. In NGC 3310 the regions 
are the central 20" and 60". The last line for each galaxy is the total flux of the 
whole galaxy. Note that the HI maps cover a much larger area than that of CO 
and SCUBA 850 /im , thus the GD mass ratio is an upper limit. 
b The 12 CO(1-0) and SCUBA 850 /im fluxes have been convolved to match the 
resolution of the HI data. The 850 ^m fluxes have been corrected for 12 CO(3-2) 
contamination, but not for thermal radio emission because its spatial distribution 
is uncertain. 

c The molecular gas mass derived using the X factor in column 6. 

d Our best estimate of the X factor in the corresponding region, in the unit of 

cm~ 2 [K kms -1 ] -1 . 

e The dust mass derived based on our best estimate of T^ ust with (3 = 2 (NGC 157) 
or an enhanced VSG model (NGC 3310) in the corresponding region (see text for 
details). 

f The 850 /jm flux here is mostly in the 60" region with SCUBA detection, but 
the HI extended to several arcminutes thus the Sg5o is a lower limit and the GD 
mass ratio is an upper limit. 



Table 10: Fluxes used for the calculation of the interstellar radiation field 





NGC3310 






NGC 157 




A 


Flux density 


Ref 


A 


Flux density 


Ref. 


[A] 


[mJy] 




[A] 


[mJy] 




1550 


90± 56 


(1) 


1516 


26±2.4 


(2) 


1910 


74 ±33 


(1) 








2267 


80± 2 


(2) 


2267 


52±2.3 


(2) 


2460 


79 ± 26 


(1) 








3320 


125±14 


(1) 








3660 (U) 


115 ±9 


(3) 


3660 (U) 


114±8 


(3) 


4400 (B) 


178±14 


(3) 


4400 (B) 


247±20 


(3) 


5530 (V) 


202 ±16 


(3) 


5530 (V) 


332± 29 


(3) 


6373 (R) 


200 ± 7 


(4) 


6400(R) 


302 ± 8 


(5) 



Notes. The fluxes are corrected for Galactic and internal extinction by extrap- 
olating the extinction in the V-band of 0.16 mag (for NGC 3310) and 0.31 (for 
NGC 157) from NED to the other wavelengths using the Draine (2003) extinction 
curve. The references for the flux densities are: (1) Code & Welch (1982), 1982, 
(2) Our own value derived from data from the GALEX archive, (3) de Vaucouleur 
et al. (1991) (RC3), (4) P. A. James, et al. ( 2004) (5) Kassin, de Jong, & Pogge 
(2006). 



